Self-Organized Criticality Driven by Deterministic Rules 
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We have investigated the essential ingredients allowing a system to show Self Organized Critical- 
ity (SOC) in its collective behavior. Using the Bak-Sneppen model of biological evolution as our 
paradigm, we show that the random microscopic rules of update can be effectively substituted with 
a chaotic map without changing the universality class. Using periodic maps SOC is preserved, but 
in a different universality class, as long as the spectrum of frequencies is broad enough. 
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Complex extended systems showing critical behavior, 
a lack of scale in their features, appear to be widespread 
in nature, being as diverse as earthquakes H, creep phe- 
nomena 0, material fracturing &-@l, fluid displacement 
in porous media WM, interface growth |p|,p[ , river net- 
works [IQ-OI and biological evolution ]l3|-[l5||. At vari- 
ance with equilibrium statistical mechanics, these sys- 
tems do not need any fine tuning of a parameter to be in 
a critical state. To explain this behavior, Bak, Tang and 
Wiesenfeld introduced the concept of self-organized criti- 
cality (SOC) through the simple sand-pile model |lq,0|. 

In recent years, several models with extremal dy- 
namics have been shown to exhibit SOC when noise is 
present 11^ . In this Letter we show that for this class of 
systems, noise can be replaced by either a chaotic or a 
quasi-periodic signal without destroying criticality. To il- 
lustrate this point we consider, as an example, the model 
proposed by Bak and Sneppen (BS) to describe the co- 
evolution of natural species |1^. The result that dif- 
ferent, even deterministic, microscopic rules can induce 
SOC in the collective behavior of a population, points to 
a greater relevance of SOC in nature. 

In the BS model an ecosystem is described by a one- 
dimensional lattice, every site of which is occupied by a 
species. Species with stronger mutual interactions in the 
ecosystem are arranged on nearest neighbor sites (the lat- 
tice can be interpreted as a food-chain, or as a food-web 
in more than one dimension). Each species is character- 
ized by its fitness, describing the average number of off- 
springs an individual of that species can have in the given 
environment. This definition of the fitness also accounts 
for the greater resistance to mutations of fitter species 
since mutations must propagate over a greater number 
of individuals to become a genetic trait of the species. 
Thus the species with the lowest fitness is the one that 
feels the strongest evolutionary pressure. Its fate is to 
either evolve or get extinct, and its place will be taken 
by some newcomer species in the same ecological niche. 
Therefore the fitness of the species occupying that site 
is the most likely to change in a short time. The near- 
est neighbor species will find a different environment, and 
their fitnesses will result changed too. As a result of such 
a simple dynamical rule, the system exhibits sequences of 
causally connected evolutionary events called avalanches 



|l8| . The number of avalanches N follows a power law 
distribution 



N{s) 



(1) 



where s is the size of the avalanche and r ^ 1.07 [|l9|J20[ 
is the avalanche critical exponent. This kind of behav- 
ior, which is the essence of self-organized criticality, has 
actually been observed in paleontological data [|l3[ sug- 
gesting that evolution and extinction may be episodic at 
all scales (a feature that goes under the name of punctu- 
ated equilibrium) p4| , |l5| . 

In nature, the evolution of the least fit species is due 
to genetic mutation. In the BS model, this mutation is 
realized by giving to the corresponding species a random 
fitness. As shown in jlq], each lattice site j is assigned a 
fitness, namely a random number between and 1. At 
each time step in the simulation the smallest fitness is 
found. Then the fitnesses of the minimum and of the 
two nearest neighbors are updated according to the rule 



fn+l — F{fn) 



(2) 



that assigns a new fitness fn+i at time n + 1 to the chosen 
lattice site. Indeed, in the original BS model, the function 
F is just a random function with a uniform distribution 
between and 1. The system reaches a stationary critical 
state in which the distribution of fitnesses is zero below a 
certain threshold fc ^ 0.66702 ll9,EQ] and uniform above 
it. It is also possible to define other quantities that show 
a power law behavior with their own critical exponent. 
Prominent among them are the first and all return time 
distributions of activity (a site is defined as active when 
its fitness is the minimum one), 



Pfit)^t- 
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Pa{t)^t- 



(3) 



where t/ - 1.58 and Ta - 0.42 |ll9|,|20|l. 

Different microscopic rules have been proposed to de- 
scribe how the mutation of the least fit species induces 
mutations of the nearest neighbors [El[ : changing the mi- 
croscopical dynamical rule affects the universality class of 
the model, but not its SOC property. This fact shows the 
robustness of the SOC behavior and poses the question 
of what are the minimal requirements the BS model has 
to satisfy in order to be critical. 



While in the BS model the updating rule (^ is a ran- 
dom function, one could consider, instead, deterministic 
updating. Indeed, we started by considering a determin- 
istic rule, whose statistical properties resemble those of 
a stochastic function, namely the Bernoulli map p2| 



fn+l = Grifn) = [r fn] , 



(4) 



where [/] stands for the value of / modulus 1 and r G N 
is a constant. It has been shown (see |22] and references 
therein) that this map has a uniform invariant measure 
(for any integer value of r) and that the Lyapunov ex- 
ponent A is given by A = logr. In Fig. |l| we show the 
power-law behavior of the first and all return probability 
distributions in the case r = 2. The critical exponents 
obtained coincide with those found in ||l^ for the ran- 
dom updating. Moreover, our simulations show that for 
all values of r the systems fall in the BS universality class, 
i.e. they all have the same critical exponents. 
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FIG. 1. First and all return distributions for a BS model 
with Bernoulli updating rule with r = 2. For all the simula- 
tions shown here, we used a lattice of 2^* sites and 5 x 10^ 
iterations exploiting the tree-algorithm explained in [23]. 

The stationary distribution of the fitnesses, on the 
other hand, follows a different pattern. Indeed, Fig. 
shows that the threshold for r = 2 is bigger than the one 
found for the random case. On increasing the value of r, 
the threshold moves towards the BS value (see Fig. ||). 
For non integer values of r (r > 1), SOC is still preserved 
within the BS universality class. However, in this case, 
the distribution of the generated numbers is not uniform 
and consequently it influences the distribution of the fit- 
nesses at the stationary state. 

The next step is then to consider updating rules that 
can be tuned to chaotic behavior by changing a parame- 
ter. 
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FIG. 2. Distribution of the fitnesses for r = 2, 3, 7, 10; 
the threshold for r = 2 is quite different from the usual BS 
threshold while the threshold corresponding to r = 10 is very 
close to the BS value (given by the vertical line). For all 
the simulations shown here, we used a lattice of 2^* sites and 
5 X 10® iterations. 

To that effect, we take as updating rule for the fitnesses 
the logistic or Feigenbaum map, namely 



/n+i — A/„(l — /„) . 



(5) 



The reasons for studying this rule are manifold. On the 
one hand, this map has already been considered in the 
context of biological evolution models and population dy- 
namics [24-26 and can thus provide a possible determin- 



istic interpretation of the evolution inside every ecological 
niche. Moreover, it has been shown that it describes the 
behavior of a wide variety of systems in nature [^ . On 
the other hand, it has a regime in which it is chaotic as 
well as one in which it is not, depending on whether A is 
bigger or less than the critical value Aoo "^ 3.56994 [ p2| . 
As Fig. y shows, for those values of A for which the map is 
chaotic, the system not only exhibits SOC but also stays 
in the same universality class as the original BS model. 
For A < Aoo we find that the system is not critical any 
more. Indeed below Aoo every site follows a periodic or- 
bit and to understand this loss of criticality one needs to 
investigate the case of periodic updating rules in every 
site. 

Let us then consider a model in which the choice of the 
new fitness is done according to the map 



fn+l — 



sin(arcsin{2fn 



1) 



1 



sin{u>j{n + 1) + ipj) + 1 



(6) 



where the loj, 0j's are, in principle, different for each site 
j and n is the time step. In our calculations we have 
chosen the frequencies tOj such that 



and the phases (^j to verify 

(j)j ^ (j)i. 



(7) 



(8) 



Bearing in mind that these maps are periodic, let us 
turn to Fig. where the first and all return probabilities 
are shown. The universality class changes with respect 
to the original BS model, with t/ = 1.67(1) and Tq = 
0.33(1), but the SOC behavior is preserved. 
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FIG. 3. First and all return distributions for a BS model 
with the logistic map with A = 4 as updating rule. The 
exponents are the same as for the BS model r/ = 1.67 and 
Ta — 0.33. In all the simulations shown in this figure, we used 
a lattice of 2^^ sites and 10* iterations. 

If instead of choosing all the frequencies different we 
choose ojj — uj and all phases different (j>j ^ (j>i the SOC 
is destroyed (even if the fitnesses are organized above a 
threshold) . This result sheds light on the loss of critical- 
ity for model (H) when A < Aoo 7 because there the system 
reduces to a set of oscillators, all with the same frequency 
and random initial phases. Model (g) should recover the 
SOC behavior if A is allowed to vary from site to site 
(mimicking condition (0) of model (||)). 

It is worth noticing that even though chaotic nraps 
produce series of numbers that may (statistically) re- 
semble random numbers (with the exception of the func- 
tional form of the invariant density), the behavior of the 
systems feels the details of the underlying dynamics, as 
shown by the dependence of the threshold on the param- 
eter r of the Bernoulli map. In particular, the Bernoulli 
map for r = 2 is formally equivalent to a coin toss [ p2|j28| 
(the paradigm of randomness) and still the threshold is 
different. However, at variance with the case of pure ran- 
dom noise, in the Bernoulli map the correlations decay 
exponentially instead of being a delta function. 
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FIG. 4. (a) First and all return distributions for a BS 
model with disordered periodic updating rules; the exponents 
are r/ = 1.67(1) and Ta = 0.33(1). (b) Distribution of the 
fitnesses. In all the simulations shown in this picture, we 
used a lattice of 2^^ sites and 5 x 10^ iterations. 

Then the system feels the details of the rule and while 
staying in the same universality class exhibits a differ- 
ent threshold. One can then conclude that, provided 
the statistical properties are those of a random process, 
SOC will persist. These results are complementary of 
the ones obtained in Ref. P9 30|: They showed that a 
chaotic system can be used instead of a heath bath to ob- 
tain thermalization (the chaotic system is referred to as 
a "booster" [^,^). This is analogous to what happens 
in the BS model: Noise (thermal or otherwise) can be re- 
placed by a deterministic chaotic system without signif- 
icant changes in the stationary state. However stochas- 
ticity in the updating rule is sufficient but not necessary: 
SOC persists even in the absence of chaos, for periodic 
updating rules. 

Summarizing, the results shown here indicate that the 
feature that ensures SOC in systems with extremal dy- 
namics, is not the randomness of the actual updated 
value but the fact that the choice of the site where the 
change is going to be performed (namely the minimum 
rule) is random. Moreover, as long as there is enough 
diversity among the species on the lattice, the longer the 
memory (or the internal correlation) of each member, the 
higher the threshold. Indeed, in the case of chaotic maps, 
the diversity is ensured by the random assignation of the 
initial values and as much as the chaoticity is increased 
we see that the threshold decreases. The extremal case 
being given by the BS model. 

In the case of the periodic map instead, the random ini- 
tial conditions are not enough to ensure diversity. Thus, 
in order to have SOC, we have to choose at random also 
the internal time-scale i. e. the periods. 

These results add strength to the relevance of SOC 



in physics and biology, since they allow different micro- 
scopic mechanisms to underlie its appearance as a collec- 
tive behavior. From the point of view of biological evolu- 
tion, this result could also account for less wild variations 
of the fitnesses. 

At this stage, several questions arise. What is the 
mechanism that allows the system to recognize the under- 
lying dynamics? Moreover, when the microscopic rule is 
disordered periodic, does the distribution of the frequen- 
cies influence the universality class of the model? 

These questions are, we think, of utmost importance 
in order to understand the widespread appearance of self 
organized criticality in nature. 
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Note: After finishing this work we became aware of 
reference |31[| where a similar line of reasoning is pursued 
for the effects of disorder on a population of integrate- 
and-fire oscillators. 
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